3. Interacting Fields 



The free field tfieories tfiat we've discussed so far are very special: we can determine 
their spectrum, but nothing interesting then happens. They have particle excitations, 
but these particles don't interact with each other. 

Here we'll start to examine more complicated theories that include interaction terms. 
These will take the form of higher order terms in the Lagrangian. We'll start by asking 
what kind of small perturbations we can add to the theory. For example, consider the 
Lagrangian for a real scalar field, 

£ = ia,0a'^0-^mV-E^'^" (3-1) 

n>3 

The coefficients A„ are called coupling constants. What restrictions do we have on A„ 
to ensure that the additional terms arc small perturbations? You might think that wc 
need simply make "A„ ^ 1". But this isn't quite right. To sec why this is the case, let's 
do some dimensional analysis. Firstly, note that the action has dimensions of angular 
momentum or, equivalently, the same dimensions as h. Since we've set h — 1, using 
the convention described in the introduction, we have [S] — 0. With S — J d'^xC, and 
[d'^x] — —4, the Lagrangian density must therefore have 

[C] = 4 (3.2) 

What does this mean for the Lagrangian (3.1)? Since [d^] — 1, we can read off the 
mass dimensions of all the factors to find, 

[0] = 1 , [m] = l , [A„] = 4-n (3.3) 

So now we see why we can't simply say we need A„ <^ 1, because this statement only 
makes sense for dimcnsionless quantities. The various terms, parameterized by A„, fall 
into three different categories 

• [A3] = 1: For this term, the dimcnsionless parameter is X3/E, where E has 
dimensions of mass. Typically in quantum field theory, E is the energy scale of 
the process of interest. This means that A3 073! is a small perturbation at high 
energies E ^ X^, but a large perturbation at low energies <^ A3. Terms that 
we add to the Lagrangian with this behavior are called relevant because they're 
most relevant at low energies (which, after all, is where most of the physics we see 
lies). In a relativistic theory, E > m, so we can always make this perturbation 
small by taking A3 <^ m. 
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• [A4] = 0: this term is small if A4 <^ 1. Such perturbations are called marginal. 

• [A„] < for n > 5: The dimensionless parameter is (A„£'"~^), which is small at 
low-energies and large at high energies. Such perturbations are called irrelevant. 

As you'll see later, it is typically impossible to avoid high energy processes in quantum 
field theory. (We've already seen a glimpse of this in computing the vacuum energy). 
This means that we might expect problems with irrelevant operators. Indeed, these lead 
to "non-renormalizable" field theories in which one cannot make sense of the infinities 
at arbitrarily high energies. This doesn't necessarily mean that the theory is useless; 
just that it is incomplete at some energy scale. 

Let me note however that the naive assignment of relevant, marginal and irrelevant 
is not always fixed in stone: quantum corrections can sometimes change the character 
of an operator. 

An Important Aside: Why QFT is Simple 

Typically in a quantum field theory, only the relevant and marginal couplings are 
important. This is basically because, as we've seen above, the irrelevant couplings 
become small at low-energies. This is a huge help: of the infinite number of interaction 
terms that we could write down, only a handful are actually needed (just two in the 
case of the real scalar field described above) . 

Let's look at this a little more. Suppose that we some day discover the true su- 
perduper "theory of everything unimportant" that describes the world at very high 
energy scales, say the GUT scale, or the Planck scale. Whatever this scale is, let's call 
it A. It is an energy scale, so [A] = 1. Now we want to understand the laws of physics 
down at our puny energy scale E <^ A. Let's further suppose that down at the energy 
scale E, the laws of physics are described by a real scalar field. (They're not of course: 
they're described by non-Abelian gauge fields and fermions, but the same argument 
applies in that case so bear with me). This scalar field will have some complicated 
interaction terms (3.1), where the precise form is dictated by all the stuff that's going 
on in the high energy supcrdupcr theory. What arc these interactions? Well, we could 
write our dimensionful coupling constants A„ in terms of dimensionless couplings Qn, 
multiplied by a suitable power of the relevant scale A, 

^n = ^, (3.4) 

The exact values of dimensionless couplings Qn depend on the details of the high-energy 
supcrduper theory, but typically one expects them to be of order 1: ~ This 
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means that for experiments at small energies E <^ A, the interaction terms of the 
form 0" with n > 4 will be suppressed by powers of [E/A)"'~^. This is usually a 
suppression by many orders of magnitude, (e.g for the energies E explored at the 
LHC, E/Mpi ~ lO"-*^^). It is this simple argument, based on dimensional analysis, that 
ensures that we need only focus on the first few terms in the interaction: those which 
are relevant and marginal. It also means that if we only have access to low-energy 
experiments (which we do!), it's going to be very difficult to figure out the high energy 
theory (which it is!), because its effects are highly diluted except for the relevant and 
marginal interactions. The discussion given above is a poor man's version of the ideas 
of effective field theory and Wilson's renormalization group, about which you can learn 
more in the "Statistical Field Theory" course. 

Examples of Weeikly Coupled Theories 

In this course we'll study only weakly coupled field theories i.e. ones that can truly be 
considered as small perturbations of the free field theory at all energies. In this section, 
we'll look at two types of interactions 

1) (j)'^ theory: 

£ = \d,ct>d^cj> - ^m'<p' - ^<P' (3.5) 

with A ^ 1. We can get a hint for what the effects of this extra term will be. Expanding 
out 4>* in terms of a- and at, we see a sum of interactions that look like 

atatatat and alalalaf; etc. (3.6) 
p p p p p p p p ^ ' 

These will create and destroy particles. This suggests that the 0^ Lagrangian describes 
a theory in which particle number is not conserved. Indeed, we could check that the 
number operator N now satisfies [H, N] ^ 0. 

2) Scalair Yukawa Theory 

L = d^ip*d^ip + ld^,(j)d^(j) - M^ip^ip - ^m^^ - ^^^^ (3.7) 

with g <^ M,m. This theory couples a complex scalar to a real scalar 0. While 
the individual particle numbers of ip and (j) are no longer conserved, we do still have 
a symmetry rotating the phase of ip, ensuring the existence of the charge Q defined 
in (2.75) such that [Q, H] = 0. This means that the number of ip particles minus the 
number of ip anti-particles is conserved. It is common practice to denote the anti- 
particle as ■0. 
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The scalar Yukawa theory has a shghtly worrying aspect: the potential has a stable 
local minimum at = ^ = 0, but is unbounded below for large enough —gcf). This 
means we shouldn't try to push this theory too far. 

A Comment on Strongly Coupled Field Theories 

In this course we restrict attention to weakly coupled field theories where we can use 
perturbative techniques. The study of strongly coupled field theories is much more 
difficult, and one of the major research areas in theoretical physics. For example, some 
of the amazing things that can happen include 

• Chcirge Fractionalization: Although electrons have electric charge 1, under 
the right conditions the elementary excitations in a solid have fractional charge 

(where N e 2Z + 1). For example, this occurs in the fractional quantum 
Hall effect. 

• Confinement: The elementary excitations of quantum chromodynamics (QCD) 
are quarks. But they never appear on their own, only in groups of three (in a 
baryon) or with an anti-quark (in a meson). They are confined. 

• Emergent Space: There are field theories in four dimensions which at strong 
coupling become quantum gravity theories in ten dimensions! The strong cou- 
pling effects cause the excitations to act as if they're gravitons moving in higher 
dimensions. This is quite extraordinary and still poorly understood. It's called 
the AdS/CFT correspondence. 

3.1 The Interaction Picture 

There's a useful viewpoint in quantum mechanics to describe situations where we have 
small perturbations to a well-understood Hamiltonian. Let's return to the familiar 
ground of quantum mechanics with a finite number of degrees of freedom for a moment. 
In the Schrodinger picture, the states evolve as 



while the operators Os are independent of time. 

In contrast, in the Heisenberg picture the states are fixed and the operators change 
in time 



L 



dt 



s 



(3.8) 



-iHt 



1^) 



H ~ 




(3.9) 
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The interaction picture is a hybrid of the two. We spht the Hamiltonian up as 

H^Ho + H,^t (3.10) 

The time dependence of operators is governed by Hq, while the time dependence of 
states is governed by i^int- Although the split into Hq and Hi^t is arbitrary, it's useful 
when Hq is soluble (for example, when Hq is the Hamiltonian for a free field theory). 
The states and operators in the interaction picture will be denoted by a subscript / 
and are given by. 



s 



Oj{t) = e^-^^^Ose-^^"* (3.11) 

This last equation also applies to i^int, which is time dependent. The interaction 
Hamiltonian in the interaction picture is. 

Hi = {H,^t)i = e'''"\H,^,)s e"^^"* (3.12) 

The Schrodinger equation for states in the interaction picture can be derived starting 
from the Schrodinger picture 



'"^^-^-Hs^s i| (e-^°*|V^),) = (i/o + i^int)5e-^^°*|^), 

^ ^ e^''°*(^int)5 e-^°* (3.13) 



So we learn that 



- Hi{t) (3.14) 



3.1.1 Dyson's Formula 



"Well, Birmingham has much the best theoretical physicist to work with, 
Peierls; Bristol has much the best experimental physicist, Powell; Cam- 
bridge has some excellent architecture. You can make your choice." 

Oppenheimer's advice to Dyson on which university position to accept. 
We want to solve (3.14). Let's write the solution as 

\m)i = u{t,tQ)mQ))j (3.15) 
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where U{t, to) is a unitary time evolution operator such that U{ti,t2)U{t2, t^) = U{ti, t^) 
and U{t,t) = 1. Then the interaction picture Schrodinger equation (3.14) requires that 

= Hr{t) U (3.16) 
If Hi were a function, then we could simply solve this by 

U(t, to) = exp (^-i Hi{t') dt'^ (3.17) 

But there's a problem. Our Hamiltonian Hj is an operator, and we have ordering 
issues. Let's see why this causes trouble. The exponential of an operator is defined in 
terms of the expansion, 

exp (^-i Hi{t') dt'^ ^l-i Hi{t') dt' + (^j^ Hi{t') dt'^ +... (3.18) 

But when we try to differentiate this with respect to t, we find that the quadratic term 
gives us 

Now the second term here looks good, since it will give part of the Hi[t)U that we 
need on the right-hand side of (3.16). But the first term is no good since the Hi{t) 
sits the wrong side of the integral term, and we can't commute it through because 
[Hj{t'),Hi{t)] 7^ when t' ^ t. So what's the way around this? 

Claim: The solution to (3.16) is given by Dyson's Formula. (Essentially first figured 
out by Dirac, although the compact notation is due to Dyson). 

C/(t, to) = T exp ^-i Hiit') dt'^ (3.20) 

where T stands for time ordering where operators evaluated at later times are placed 
to the left 

T(O,(t0(^2(t.)) = l^^l'^J^i'^J '''''' (3.21) 
\02{t2)0^{t^) t2>t. 

Expanding out the expression (3.20), we now have 

U{t, to)^l-i f dt' Hi{t') + f dt' f dt" Hi{t")Hi{t') 

+ f dt' f dt" Hi{t')Hi{t") 



+ 
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Actually these last two terms double up since 

pt ft ft" 

dt' / dt" Hi{t")Hi{t') = / dt" / dt' Hi{t")Hi{t') 

Jt' Jto Jto 



t dt' t dt" Hi{t')Hi{t") (3.22) 

J to J to 



where the range of integration in the first expression is over t" > t', while in the second 
expression it is t' < t" which is, of course, the same thing. The final expression is the 
same as the second expression by a simple relabelling. This means that we can write 



U{t,tQ) ^l-i I dt'Hj(t') + (-if I dt' I dt" Hi(t')Hj(t" 

''to J to J to 



) + ... (3.23) 



Proof: The proof of Dyson's formula is simpler than explaining what all the nota- 
tion means! Firstly observe that under the T sign, all operators commute (since their 
order is already fixed by the T sign) . Thus 

T exp (^-i dt'Hiit')^ = T Hi{t) ex.p (^-i dt' Hi{t')^ 

= Hi{t) T exp ^ dt' Hi{t')^ (3.24) 

since t, being the upper limit of the integral, is the latest time so Hi{t) can be pulled 
out to the left. □ 

Before moving on, I should confess that Dyson's formula is rather formal. It is 
typically very hard to compute time ordered exponentials in practice. The power of 
the formula comes from the expansion which is valid when Hj is small and is very easily 
computed. 

3.2 A First Look at Scattering 

Let us now apply the interaction picture to field theory, starting with the interaction 
Hamiltonian for our scalar Yukawa theory, 

Hi^t ^gjd^x V^t^0 (3.25) 

Unlike the free theories discussed in Section 2, this interaction doesn't conserve particle 
number, allowing particles of one type to morph into others. To see why this is, we use 
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the interaction picture and follow the evolution of the state: {ipit)) = f/(t, to) l'^(^o)); 
where U{t,to) is given by Dyson's formula (3.20) which is an expansion in powers of 
i^int- But i^int contains creation and annihilation operators for each type of particle. 
In particular, 

• ~ a + a^: This operator can create or destroy (p particles. Let's call them 
mesons. 

• ijj ~ 6 + c^: This operator can destroy particles through b, and create anti- 
particles through c^. Let's call these particles nucleons. Of course, in reality 
nucleons are spin 1/2 particles, and don't arise from the quantization of a scalar 
field. But we'll treat our scalar Yukawa theory as a toy model for nucleons 
interacting with mesons. 

• ip^ b"^ -\- c: This operator can create nucleons through b^, and destroy anti- 
nucleons through c. 

Importantly, Q = — remains conserved in the presence of i/int- At first order in 
perturbation theory, we find terms in Hi^t like c^b^a. This kills a meson, producing a 
nucleon-anti-nucleon pair. It will contribute to meson decay 4> 4>'4^- 

At second order in perturbation theory, we'll have more complicated terms in (ifint)^, 
for example {c^b'^a){cba^). This term will give contributions to scattering processes 
ipip — > ^ ijjilj. The rest of this section is devoted to computing the quantum ampli- 
tudes for these processes to occur. 

To calculate amplitudes we make an important, and slightly dodgy, assumption: 

Initial and final states are eigenstates of the free theory 

This means that we take the initial state \i) at t — >• — oo, and the final state |/) at 
t — > -|-oo, to be eigenstates of the free Hamiltonian Hq. At some level, this sounds 
plausible: at i — > — oo, the particles in a scattering process are far separated and don't 
feel the effects of each other. Furthermore, we intuitively expect these states to be 
eigenstates of the individual number operators N , which commute with ifo, but not 
i^int- As the particles approach each other, they interact briefiy, before departing again, 
each going on its own merry way. The amplitude to go from \i) to |/) is 

hm (/|C/(t+,t_)|i) = (/|5|i) (3.26) 

t±— >±oo 

where the unitary operator S is known as the S-matrix. (S is for scattering). There 
are a number of reasons why the assumption of non-interacting initial and final states 
is shaky: 
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• Obviously we can't cope with bound states. For example, this formalism can't 
describe the scattering of an electron and proton which collide, bind, and leave 
as a Hydrogen atom. It's possible to circumvent this objection since it turns out 
that bound states show up as poles in the S-matrix. 

• More importantly, a single particle, a long way from its neighbors, is never alone 
in field theory. This is true even in classical electrodynamics, where the electron 
sources the electromagnetic field from which it can never escape. In quantum 
electrodynamics (QED), a related fact is that there is a cloud of virtual photons 
surrounding the electron. This line of thought gets us into the issues of renormal- 
ization — more on this next term in the "AQFT" course. Nevertheless, motivated 
by this problem, after developing scattering theory using the assumption of non- 
interacting asymptotic states, we'll mention a better way. 

3.2.1 An Example: Meson Decay 

Consider the relativistically normalized initial and final states, 

|i) = y/2E^al |0) 

= V4£p^feJ.4 |0) (3.27) 

The initial state contains a single meson of momentum p; the final state contains a 
nucleon-anti-nucleon pair of momentum qi and q2. We may compute the amplitude for 
the decay of a meson to a nucleon-anti-nucleon pair. To leading order in g, it is 

= -ig{f\ J d^x^\x)^{x)(t>{x) |i) (3.28) 

Let's go slowly. We first expand out ~ a + a''" using (2.84). (Remember that the (j) in 
this formula is in the interaction picture, which is the same as the Heisenberg picture 
of the free theory). The a piece will turn \i) into something proportional to |0), while 
the piece will turn \i) into a two meson state. But the two meson state will have 
zero overlap with (/|, and there's nothing in the -0 and operators that lie between 
them to change this fact. So we have 

{f\S\i) = -ig{f\ j d'xi;\x)i;{x) I (03^^«fc-4^"''''|O) 

^-igifl j d'^x^p\x)^p{x)e-'P-''\Q) (3.29) 

where, in the second line, we've commuted past at, picking up a 6^^\p- k) delta- 
function which kills the d^k integral. We now similarly expand out ~ 6 -|- and 
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ip^ + c. To get non-zero overlap with (/|, only the and contribute, for they 
create the nucleon and anti-nucleon from |0). We then have 



= -ig{2n)U^^\q, + q2-p) (3.30) 
and so we get our first quantum field theory amplitude. 

Notice that the 5-function puts constraints on the possible decays. In particular, the 
decay only happens at all if m > 2M. To see this, we may always boost ourselves 
to a reference frame where the meson is stationary, so p = (m, 0,0, 0). Then the 
delta function imposes momentum conservation, telling us that qi — —q2 and m — 

Later you will learn how to turn this quantum amplitude into something more phys- 
ical, namely the lifetime of the meson. The reason this is a httle tricky is that we must 
square the amplitude to get the probability for decay, which means we get the square 
of a 5-function. We'll explain how to deal with this in Section 3.6 below, and again in 
next term's "Standard Model" course. 

3.3 Wick's Theorem 

Prom Dyson's formula, we want to compute quantities like (/| T {Hi{xi) . . . Hj{xn)} 
where \i) and |/) are eigenstates of the free theory. The ordering of the operators is fixed 
by T, time ordering. However, since the i?/'s contain certain creation and annihilation 
operators, our life will be much simpler if we can start to move all annihilation operators 
to the right where they can start killing things in Recall that this is the definition 
of normal ordering. Wick's theorem tells us how to go from time ordered products to 
normal ordered products. 

3.3.1 An Example: Recovering the Propagator 

Let's start simple. Consider a real scalar field which we decompose in the Heisenberg 
picture as 

(P{x) = (p-^{x) + (p-{x) (3.31) 



where 



p 

p 



(2ir)3 
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where the dz signs on 0^ make Uttle sense, but apparently you have Pauh and Heisenberg 
to blame. (They come about because 0+ ~ e~*^*, which is sometimes called the positive 
frequency piece, while 0~ ~ Q+t^t ^^^q negative frequency piece). Then choosing 
x° > ?/°, we have 

T(l>{x)cj>{y)^(l>{x)cf>{y) 

= (0+(x) + 0-(x))(0+(y) + r{y)) (3.33) 

= 0+(x)0+(y) + 0-(x)0+(y) + 0-(y)0+(x) + [0+(x), riv)] + r{x)r{y) 

where the last line is normal ordered, and for our troubles we have picked up a factor 
of D{x — y) = [(j)^{x), (t)~{y)] which is the propagator we met in (2.90). So for x^ > y° 
we have 

T (t>{x)ct>{y) =: 0(x)0(y) :+D{x-y) (3.34) 
Meanwhile, for 7/° > we may repeat the calculation to find 

T (t>{x)cf>{y) =: 0(x)0(y) :+£>(?/- x) (3.35) 
So putting this together, we have the final expression 

T ct>{x)ct>{y) =: (t>{x)ct>{y) : + (x - y) (3.36) 

where Ai7(a; — y) is the Feynman propagator defined in (2.93), for which we have the 
integral representation 

/j4l -ik-ix-y) 
[2n)* — + le 

Let me reiterate a comment from Section 2: although T (j){x)(f){y) and : 0(a;)0(|/) : are 
both operators, the difference between them is a c- number function, Af{x — y). 

Definition: We define the contraction of a pair of fields in a string of operators 
. . . (t>{xi) . . . 4>{x2) ... to mean replacing those operators with the Feynman propaga- 
tor, leaving all other operators untouched. We use the notation. 



...0(xi)...0(x2)... (3.38) 
to denote contraction. So, for example, 

'^{x)^ = Ap{x - y) (3.39) 
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A similar discussion holds for complex scalar fields. We have 

T'il;{x)iP\y) =: i^ix)ip''{y) : +Af{x - y) (3.40) 
prompting us to define the contraction 

= Ai.(x -y) and = V(2;)V(z/)' = (3.41) 

3.3.2 Wick's Theorem 

For any collection of fields 0i = (p{xi), 02 = 0(2^2), etc, we have 

T(0i . . . 0„) =: 01 . . . 0„ : + : all possible contractions : (3.42) 

To see what the last part of this equation means, let's look at an example. For n — A, 
the equation reads 

T(0i020304) = : 01020304 : +0102 : 0304 : +0103 '■ 0204 '■ + four similar terms 

+ 0102 0304 + 0103 0204 + 0104 0203 (3.43) 

Proof: The proof of Wick's theorem proceeds by induction and a little thought. It's 
true for n = 2. Suppose it's true for 02 • • • 0„ and now add 0i. We'll take x^ > x^ 
for all A; = 2, . . . , n. Then we can pull 0i out to the left of the time ordered product, 
writing 

T (0102 . . . 0n) = (0^ + 0r) (: 02 ■ ■ ■ 0n : + : coutractious :) (3.44) 

The (pi term stays where it is since it is already normal ordered. But in order to write 
the right-hand side as a normal ordered product, the (pf term has to make its way 
past the crowd of 0^ operators. Each time it moves past 0^, we pick up a factor of 

0i0fc = Af(xi — Xk) from the commutator. (Try it!) □ 

3.3.3 An Example: Nucleon Scattering 

Let's look at ipip ipip scattering. We have the initial and final states 

\i) = y/2Ei,y/2Ei,bl bl |0) = |pi,p2) 

I/) = ^/2E^^^/2E^,bl,bl, |0) = \p[,p^) (3.45) 

We can then look at the expansion of {f\S\i). In fact, wc really want to calculate 
{f\S— since we're not interested in situations where no scattering occurs. At 
order we have the term 

^ d^x,d''x2 T (V'^(xi)V'(xi)0(xi)V'^(x2)V(2:2)0(2:2)) (3.46) 
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Now, using Wick's theorem we see there is a piece in the string of operators which looks 
hke 



: '4'\xi)ij{xi)ij\x2)i^{x2) : (p{xi)(f){x2) (3.47) 

which will contribute to the scattering because the two ip fields annihilate the ip parti- 
cles, while the two ip'' fields create ip particles. Any other way of ordering the ip and 
ip'^ fields will give zero contribution. This means that we have 

{p'l,P2\ ■ 1p\xi)'lp{xi)'lp''{x2)'lp{x2) : bl,P2) 

= {p'l,P2\^\xi)ij\x2) |0) {0\lp{Xi)l/j{x2) \PUP2) 

= ^e*'^'i'^'^''"'^2'^2 _j_ ^ip[-X2+ip'2-Xl^ ^^-ipi-Xl-ip2-X2 _j_ ^-ipi-X2-ip2-xi^ 

— Qixi-{p[-pi)+iX2-ip'2-P2) _|_ gia;i-(p^-pi)+ix2-(p'i-p2) _|_ ^ (3.48) 

where, in going to the third line, we've used the fact that for relativistically normalized 
states, 

(0|V'(2:) |p) = e-^f^ (3.49) 
Now let's insert this into (3.46), to get the expression for (/| S \i) at order g'^, 

(—jn\^ r r d'^h jpik-(xi-X2) 

^ / ^^.^^'^ + + « X.,] / A_ ___ (3.50, 

where the expression in square brackets is (3.48), while the final integral is the 
propagator which comes from the contraction in (3.47). Now the {xi <-> ^2) terms 
double up with the others to cancel the factor of 1/2 out front. Meanwhile, the Xi and 
X2 integrals give delta-functions. We're left with the expression 

+ •S''" W - Pi + k) W - K - 4)] (3.51) 

Finally, we can trivially do the d^k integral using the delta-functions to get 



1 1 

+ 



(pi ~ p'lY ~ + (pi ~ P2Y ~ + 



(27r)^(5(^)(pi+P2-pi-p^) 



In fact, for this process we may drop the +ie terms since the denominator is never 
zero. To see this, we can go to the center of mass frame, where pi — —p2 and, by 
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momentum conservation, = This ensures that the 4-momentum of the meson 
is A; = (0,p — p '), so A;^ < 0. We therefore have the end result, 



1 

+ 



{27i)U^'\p, + P2 - p[ - p',) (3.52) 



We will see another, much simpler way to reproduce this result shortly using Feynman 
diagrams. This will also shed light on the physical interpretation. 

This calculation is also relevant for other scattering processes, such as ^bip — > 
t/jt/j — >• ipi/j. Each of these comes from the term (3.48) in Wick's theorem. However, we 
will never find a term that contributes to scattering ipip ipip, for this would violate 
the conservation of Q charge. 

Another Example: Meson-Nucleon Scattering 

If we want to compute ilxp ^ scattering at order g'^, we would need to pick out the 
term 

: {xi)(i){xi)^{x2)(i){x2) : t/'(a;i)t/^^(a;2) (3.53) 

and a similar term with ip and exchanged. Once more, this term also contributes to 
similar scattering processes, including T/i^ -^xf) and 00 — > 

3.4 Feynman Diagrams 

"Like the silicon chips of more recent years, the Feynman diagram was 

bringing computation to the masses." , ,. „ , . 

Julian bcnwmger 

As the above example demonstrates, to actually compute scattering amplitudes using 
Wick's theorem is rather tedious. There's a much better way. It requires drawing pretty 
pictures. These pictures represent the expansion of (/| S \i) and we will learn how to 
associate numbers (or at least integrals) to them. These pictures are called Feynman 
diagrams. 

The object that we really want to compute is (/| S—1 \i), since we're not interested in 
processes where no scattering occurs. The various terms in the perturbative expansion 
can be represented pictorially as follows 

• Draw an external line for each particle in the initial state \i) and each particle 
in the final state |/). We'll choose dotted lines for mesons, and solid lines for 
nucleons. Assign a directed momentum p to each line. Further, add an arrow to 
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solid lines to denote its charge; we'll choose an incoming (outgoing) arrow in the 
initial state for ip (^). We choose the reverse convention for the final state, where 
an outgoing arrow denotes ip. 



• Join the external lines together with trivalent vertices ^ 




Each such diagram you can draw is in 1-1 correspondence with the terms in the 
expansion oi {f\ S — 

3.4.1 Feynman Rules 

To each diagram we associate a number, using the Feynman rules 

• Add a momentum k to each internal line 

• To each vertex, write down a factor of 



where ki is the sum of all momenta flowing into the vertex. 

• For each internal dotted line, corresponding to a particle with momentum k, 
we write down a factor of 



We include the same factor for solid internal ^jJ lines, with m replaced by the 
nucleon mass M. 




(3.54) 




(3.55) 
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3.5 Examples of Scattering Amplitudes 

Let's apply the Feynman rules to compute the amplitudes for various processes. We 
start with something familiar: 

Nucleon Scattering Revisited 

Let's look at how this works for the ipip — > ipip scattering at order g^. We can write 
down the two simplest diagrams contributing to this process. They are shown in Figure 
9. 




Figure 9: The two lowest order Feynman diagrams for nucleon scattering. 
Applying the Feynman rules to these diagrams, we get 



+ 



(27rr5(^)(pi + P2-p'i-p'2) (3.56) 



which agrees with the calculation (3.51) that we performed earlier. There is a nice 
physical interpretation of these diagrams. We talk, rather loosely, of the nucleons 
exchanging a meson which, in the first diagram, has momentum k — (pi—Pi) = (P2— ^2)- 
This meson doesn't satisfy the usual energy dispersion relation, because k'^ ^rr?: the 
meson is called a virtual particle and is said to be off-shell (or, sometimes, off mass- 
shell). Heuristically, it can't live long enough for its energy to be measured to great 
accuracy. In contrast, the momentum on the external, nucleon legs all satisfy = M^, 
the mass of the nucleon. They are on-shell. One final note: the addition of the two 
diagrams above ensures that the particles satisfy Bose statistics. 

There are also more complicated diagrams which will contribute to the scattering 
process at higher orders. For example, we have the two diagrams shown in Figures 
10 and 11, and similar diagrams with p[ and P2 exchanged. Using the Feynman rules, 
each of these diagrams translates into an integral that we will not attempt to calculate 
here. And so we go on, with increasingly complicated diagrams, all appearing at higher 
order in the coupling constant g. 
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I 






Figure 10: A contribution at 0{g'^ 



Figure 11: A contribution at 0{g^ 



Amplitudes 

Our final result for the nucleon scattering amplitude {f\S—l\i) at order g"^ was 



The ^-function follows from the conservation of 4-momcntum which, in turn, follows 
from spacetime translational invariancc. It is common to all S-matrix elements. We will 
define the amplitude Afi by stripping off this momentum-conserving delta-function. 



{f\S-l\i)^iAfi{2n)'S^'\pF-pi) 



(3.57) 



where pi (pp) is the sum of the initial (final) 4-momenta, and the factor of i out front 
is a convention which is there to match non-relativistic quantum mechanics. We can 
now refine our Feynman rules to compute the amplitude iAfi itself: 



Draw all possible diagrams with appropriate external legs and impose 4- momentum 
conservation at each vertex. 



• Write down a factor of {—ig) at each vertex. 

• For each internal line, write down the propagator 

• Integrate over momentum k flowing through each loop J d'^k/ {2Tr)^. 

This last step deserves a short explanation. The diagrams we've computed so far have 
no loops. They are tree level diagrams. It's not hard to convince yourself that in 
tree diagrams, momentum conservation at each vertex is sufficient to determine the 
momentum flowing through each internal line. For diagrams with loops, such as those 
shown in Figures 10 and 11, this is no longer the case. 
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Figure 12: The two lowest order Feynman diagrams for nucleon to meson scattering. 



Nucleon to Meson Scattering 

Let's now look at the amplitude for a nucleon-anti-nucleon pair to annihilate into a 
pair of mesons: — > (fxf). The simplest Feynman diagrams for this process are shown 
in Figure 12 where the virtual particle in these diagrams is now the nucleon ip rather 
than the meson 0. This fact is reflected in the denominator of the amplitudes which 
are given by 



iA 



M2 



(3.58) 



L(pi-K)2-M2 (pi-p^)2 

As in (3.52), we've dropped the ie from the propagators as the denominator never 
vanishes. 



Nucleon-Anti-Nucleon Scattering 






Figure 13: The two lowest order Feynman diagrams for nucleon-anti-nucleon scattering. 

For the scattering of a nucleon and an anti-nucleon, ip'ip — > the Feynman 
diagrams arc a little different. At lowest order, they are given by the diagrams of 
Figure 13. It is a simple matter to write down the amplitude using the Feynman rules. 



ipi — — 'm? (pi + ~ + 



(3.59) 
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(nb) 



ALEPH 



Notice that the momentum dependence in the sec- 
ond term is different from that of nucleon-nucleon 
scattering (3.56), reflecting the different Feynman di- 
agram that contributes to the process. In the center 
of mass frame, pi = — p2, the denominator of the sec- 
ond term is 4(M2 + pi) - m^. If m < 2M, then this 
term never vanishes and we may drop the ie. In con- 
trast, if m > 2M, then the amplitude corresponding 
to the second diagram diverges at some value of p. , 
In this case it turns out that we may also neglect the 
ie term, although for a different reason: the meson 
is unstable when m > 2M, a result we derived in 
(3.30). When correctly treated, this instabihty adds 
a finite imaginary piece to the denominator which 

overwhelms the ie. Nonetheless, the increase in the scattering amplitude which we see 
in the second diagram when 4(M^ +P^) = is what allows us to discover new parti- 
cles: they appear as a resonance in the cross section. For example, the Figure 14 shows 
the cross-section (roughly the amplitude squared) plotted vertically for e^e" fi^p~ 
scattering from the ALEPH experiment in CERN. The horizontal axis shows the center 
of mass energy. The curve rises sharply around 91 GeV, the mass of the Z-boson. 




Figure 14: 



Meson Scattering 

For 00 — > 00, the simplest diagram we can write 
down has a single loop, and momentum conservation at 
each vertex is no longer sufficient to determine every 
momentum passing through the diagram. We choose 
to assign the single undetermined momentum k to the 
right-hand propagator. All other momenta are then de- 
termined. The amplitude corresponding to the diagram 
shown in the figure is 



HgY 



1 



(27r)4 (k^ - + ie){{k + p[)^ - + ie) 

1 




Figure 15: 



X 



{{k + p{ - pi)2 - M2 + ie){{k - p^)2 - M2 + ie) 



These integrals can be tricky. For large k, this integral goes as J d^k/k^, which is at 
least convergent as k —>■ oo. But this won't always be the case! 
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3.5.1 Mandelstam Variables 

We see that in many of the amplitudes above — in particular those that include the 
exchange of just a single particle — the same combinations of momenta are appear- 
ing frequently in the denominators. There are standard names for various sums and 
differences of momenta: they are known as Mandelstam variables. They are 

s = (Pl+P2)' = (P'l+Ps)' 

t^{pi-p'if^{p2-P2f (3.60) 

where, as in the examples above, pi and p2 are the momenta of the two initial particles, 
and p[ and P2 are the momenta of the final two particles. We can define these variables 
whether the particles involved in the scattering are the same or different. To get a feel 
for what these variables mean, let's assume all four particles are the same. We sit in 
the center of mass frame, so that the initial two particles have four-momenta 

pi = {E, 0, 0, p) and p2 = {E, 0, 0, -p) (3.61) 

The particles then scatter at some angle 9 and leave with momenta 

p'^ = (£^, O,psin^,pcos^) and P2 = (£^, 0, — psin^, — pcos^) (3.62) 

Then from the above definitions, we have that 

s = 4:E^ and t = -2p2(l - cos 6) and u = -2p2(l + cos 6) (3.63) 

The variable s measures the total center of mass energy of the coUision, while the 
variables t and u are measures of the momentum exchanged between particles. (They 
are basically equivalent, just with the outgoing particles swapped around). Now the 

amplitudes that involve exchange of a single particle can be written simply in terms of 
the Mandelstam variables. For example, for nucleon-nucleon scattering, the amplitude 
(3.56) is schematically A ^ {t — m^)~^ + {u — m^)~^. For the nucleon-anti-nucleon 
scattering, the amplitude (3.59) is ^ ~ (t — rn^)^^ + (s — 171"^)^^. We say that the 
first case involves "t-channel" and "u-channel" diagrams. Meanwhile the nucleon-anti- 
nucleon scattering is said to involve "t-channel" and "s-channel" diagrams. (The first 
diagram indeed includes a vertex that looks like the letter "T" ) . 

Note that there is a relationship between the Mandelstam variables. When all the 
masses are the same we have s + t + u = 4M^. When the masses of all 4 particles differ, 
this becomes s + t + u = M^. 
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3.5.2 The Yukawa Potential 

So far we've computed the quantum amplitudes for various scattering processes. But 
these quantities are a httle abstract. In Section 3.6 below (and again in next term's 
"Standard Model" course) we'll see how to turn amplitudes into measurable quantities 
such as cross-sections, or the lifetimes of unstable particles. Here we'll instead show 
how to translate the amplitude (3.52) for nucleon scattering into something familiar 
from Newtonian mechanics: a potential, or force, between the particles. 

Let's start by asking a simple question in classical field theory that will turn out to 
be relevant. Suppose that we have a fixed (5-function source for a real scalar field (p, 
that persists for all time. What is the profile of (f){x)7 To answer this, we must solve 
the static Klein-Gordon equation, 

-VV + mV = 5(^)(f) (3.64) 

We can solve this using the Fourier transform. 

Plugging this into (3.64) tells us that (A;^ -|- m'^)(p{k) — 1, giving us the solution 

/j3jL pik-x 
7^ ^ (3-66) 
(27r)3A;2 + ^2 

Let's now do this integral. Changing to polar coordinates, and writing k ■ x = kr cos 6*, 
we have 

1 /"^ „ k^ 2 sin kr 



1 /"^ 



/j.2 _|_ ^2 J^j. 



1 /■+°° „ A; sin kr 
dk 



(27r)2r k"^ + m? 



^ -Re 



27rr 



dk_ 
J-oo 27ri 



^2 _|_ j^2 



(3.67) 



We compute this last integral by closing the contour in the upper half plane k +ioo, 
picking up the pole at k — +im. This gives 

0(f) = J- e-""' (3.68) 

The field dies off exponentially quickly at distances 1/m, the Compton wavelength of 
the meson. 
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Now we understand the profile of the field, what does this have to do with the force 
between ip particles? We do very similar calculations to that above in electrostatics 
where a charged particle acts as a 5-function source for the gauge potential: —V^Aq — 
5^^\x), which is solved by Aq = l/Anr. The profile for Aq then acts as the potential 
energy for another charged (test) particle moving in this background. Can we give the 
same interpretation to our scalar field? In other words, is there a classical limit of the 
scalar Yukawa theory where the particles act as 5-function sources for 0, creating 
the profile (3.68)? And, if so, is this profile then felt as a static potential? The answer 
is essentially yes, at least in the limit M ^ m. But the correct way to describe the 
potential felt by the ip particles is not to talk about classical fields at all, but instead 
work directly with the quantum amplitudes. 

Our strategy is to compare the nucleon scattering amplitude (3.52) to the corre- 
sponding amplitude in non-relativistic quantum mechanics for two particles interacting 
through a potential. To make this comparison, we should first take the non-relativistic 
limit of (3.52). Let's work in the center of mass frame, with p = pi = —p2 and 
p' = Pi = —P2- The non-relativistic limit means \p\ <^ M which, by momentum 
conservation, ensures that \p'\ <^ M. In fact one can check that, for this particular 
example, this limit doesn't change the scattering amplitude (3.52): it's given by 



1 1 



{p — p'Y + m? {p + p'Y + ni? 



(3.69) 



How do we compare this to scattering in quantum mechanics? Consider two particles, 
separated by a distance r, interacting through a potential U{r). In non-relativistic 
quantum mechanics, the amplitude for the particles to scatter from momentum states 
±p into momentum states ±p' can be computed in perturbation theory, using the 
techniques described in Section 3.1. To leading order, known in this context as the 
Born approximation, the amplitude is given by 

(p'l U{r) \p) J d^rU{f)e-'^-^>^ (3.70) 

There's a relative factor of (2M)^ that arises in comparing the quantum field theory 
amplitude A to {p'\ U{r) that can be traced to the rclativistic normalization of the 
states \piiP2)- (It is also necessary to get the dimensions of the potential to work out 
correctly). Including this factor, and equating the expressions for the two amplitudes, 
we get 

dh U(r) e-'^-^>^ = — — — (3.71) 



-68- 



where we've introduced the dimensionless parameter A — g/2M. We can trivially invert 
this to find, 



But this is exactly the integral (3.66) we just did in the classical theory. We have 



This is the Yukawa potential. The force has a range 1/m, the Compton wavelength of 
the exchanged particle. The minus sign tells us that the potential is attractive. 

Notice that quantum field theory has given us an entirely new perspective on the 
nature of forces between particles. Rather than being a fundamental concept, the force 
arises from the virtual exchange of other particles, in this case the meson. In Section 6 
of these lectures, we will see how the Coulomb force arises from quantum field theory 
due to the exchange of virtual photons. 

We could repeat the calculation for nuclcon-anti-nucleon scattering. The amplitude 
from field theory is given in (3.59). The first term in this expression gives the same 
result as for nucleon-nucleon scattering with the same sign. The second term vanishes in 
the non-rclativisitic limit (it is an example of an interaction that doesn't have a simple 
Newtonian interpretation). There is no longer a factor of 1/2 in (3.70), because the 
incoming/outgoing particles are not identical, so we learn that the potential between 
a nucleon and anti-nucleon is again given by (3.73). This reveals a key feature of 
forces arising due to the exchange of scalars: they are universally attractive. Notice 
that this is different from forces due to the exchange of a spin 1 particle — such as 
electromagnetism — where the sign flips when we change the charge. However, for 
forces due to the exchange of a spin 2 particle — i.e. gravity — the force is again 
universally attractive. 

3.5.3 0^ Theory 

Let's briefly look at the Feynman rules and scattering amphtudes for the interaction 
Hamiltonian 



The theory now has a single interaction vertex, which comes with a factor of {—iX), 
while the other Feynman rules remain the same. Note that we assign (— iA) to the 




(3.72) 




(3.74) 
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vertex rather than (— ■iA/4!). To see why this is, we can look at 00 scattering, 

which has its lowest contribution at order A, with the term 

{p[,P2\ ■ (l){x)(j){x)(j){x)(/){x) : \pi,P2) (3.75) 

Any one of the fields can do the job of annihilation or creation. This gives 4! different 
contractions, which cancels the 1/4! sitting out front. 

Feynman diagrams in the 0^ theory sometimes come with extra 
combinatoric factors (typically 2 or 4) which are known as symmetry 
factors that one must take into account. For more details, see the book 
by Peskin and Schroeder. 

Using the Feynman rules, the scattering amplitude for 00 (fxf) is Figure 16: 
simply iA = — iA. Note that it doesn't depend on the angle at which 
the outgoing particles emerge: in 0'^ theory the leading order two-particle scattering 
occurs with equal probability in all directions. Translating this into a potential between 
two mesons, we have 

So scattering in 0^ theory is due to a 5-function potential. The particles don't know 
what hit them until it's over. 

3.5.4 Connected Diagrams and Amputated Diagrams 

We've seen how one can compute scattering amplitudes by writing down all Feynman 
diagrams and assigning integrals to them using the Feynman rules. In fact, there are 
a couple of caveats about what Feynman diagrams you should write down. Both of 
these caveats are related to the assumption we made earlier that "initial and final states 
are eigenstates of the free theory" which, as we mentioned at the time, is not strictly 
accurate. The two caveats which go some way towards amehorating the problem are 
the following 

• We consider only connected Feynman diagrams, where every part of the diagram 
is connected to at least one external line. As we shall see shortly, this will be 
related to the fact that the vacuum |0) of the free theory is not the true vacuum 
1^2) of the interacting theory. An example of a diagram that is not connected is 
shown in Figure 17. 




- 70 - 



• We do not consider diagrams with loops on external lines, for example the diagram 
shown in the Figure 18. We will not explain how to take these into account in this 
course, but you will discuss them next term. They are related to the fact that the 
one-particle states of the free theory are not the same as the one-particle states of 
the interacting theory. In particular, correctly dealing with these diagrams will 
account for the fact that particles in interacting quantum field theories are never 
alone, but surrounded by a cloud of virtual particles. We will refer to diagrams 
in which all loops on external legs have been cut-off as "amputated" . 




Figure 17: A disconnected diagram. Figure 18: An un-amputated diagram 



3.6 What We Measure: Cross Sections and Decay Rates 

So far we've learnt to compute the quantum amplitudes for particles decaying or scat- 
tering. As usual in quantum theory, the probabilities for things to happen are the 
(modulus) square of the quantum amplitudes. In this section we will compute these 
probabilities, known as decay rates and cross sections. One small subtlety here is that 
the S-matrix elements (/| S' — 1 |i) all come with a factor of {2n)^S^^\pF — Pi), so we 
end up with the square of a delta-function. As we will now see, this comes from the 
fact that we're working in an infinite space. 



3.6.1 Fermi's Golden Rule 

Let's start with something familiar and recall how to derive Fermi's golden rule from 
Dyson's formula. For two energy eigenstates \m) and \n), with Em 7^ En, we have to 
leading order in the interaction, 

(m| U{t) \n) = -i (m| [ dtHi(t) \n) 

Jo 

t 



—i {m\ Hiat\n) I dt' e 
Jo 



J,ujt _ 1 

- (m| i/int \n) (3.77) 
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where uj = Em — En- This gives us the probabihty for the transition from \n) to |m) in 
time t, as 



P„(t) = I (m| Uit) \n)\^ = 2| (m| i/i^t \n) \' 

The function in brackets is plotted in Figure 19 for fixed t. 
We see that in time t, most transitions happen in a region 
between energy eigenstates separated by AE = 27i/t. As 
t — > oo, the function in the figure starts to approach a delta- 
function. To find the normalization, we can calculate 



1 — cos tut 



(3.78) 



duj 

oo 

1 — COS ut 



1 — COS cut 



TTt 




Trt6{u!) as t oo 



Figure 19: 



Consider now a transition to a cluster of states with density 
p{E). In the limit t — »• oo, we get the transition probability 



P. 



dEm p{Em) 2| (m| ifint \n) \ ^ 
2tX I (m| i/int 1^) |^p(-En)t 



1 — COS cut 



(3.79) 



which gives a constant probability for the transition per unit time for states around 
the same energy En ~ E^ = E. 



P. 



27r| {m\H,M ? p{E) 



(3.80) 



This is Fermi's Golden Rule. 



In the above derivation, we were fairly careful with taking the limit as t oo. 
Suppose we were a little sloppier, and first chose to compute the amplitude for the 
state \n) at t ^ — oo to transition to the state \m) at t ^ +oo. Then we get 

/■t=+oo 

-i (m| / Hi{t) \n) = -i {m\ Hm \n) 2tt5{uj) (3.81) 

Now when squaring the amplitude to get the probability, we run into the problem of 
the square of the delta-function: P„^m = | (^| -f^int |^) \'^{27t)'^S{ljY. Tracking through 
the previous computations, we realize that the extra infinity is coming because Pm^n 
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is the probability for the transition to happen in infinite time t ^ oo. We can write 
the delta-functions as 

(27r)^5(a;)^ = (27r)%) T (3.82) 

where T is shorthand for t ^ oo (we used a very similar trick when looking at the 
vacuum energy in (2.25)). We now divide out by this power of T to get the transition 
probability per unit time, 

P„ = 27r| (m| i/int \n) I" 5{cu) (3.83) 

which, after integrating over the density of final states, gives us back Fermi's Golden 
rule. The reason that we've stressed this point is because, in our field theory calcula- 
tions, we've computed the amplitudes in the same way as (3.81), and the square of the 
5*^^^ -functions will just be re-interpreted as spacetime volume factors. 

3.6.2 Decay Rates 

Let's now look at the probability for a single particle \i) of momentum p/ (I=initial) 
to decay into some number of particles |/) with momentum pi and total momentum 
Pf = ^iPi- This is given by 

Our states obey the relativistic normalization formula (2.65), 

(i| i) = (27r)^ 2E^^ 5^^\0) = 2Ep-V (3.85) 
where we have replaced 5*^^)(0) by the volume of 3-space. Similarly, 

(/!/)= n (3.86) 

final states 

If we place our initial particle at rest, so pV = and E^^ — m, we get the probability 
for decay 

pJ^^2n)V^>^„.,,^yT ^ ^ (3.87) 

final states 

where, as in the second derivation of Fermi's Golden Rule, we've exchanged one of the 
delta- functions for the volume of spacetime: (27r)^5*^^'*(0) = VT. The amplitudes Afi 
are, of course, exactly what we've been computing. (For example, in (3.30), we saw 
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that A = —g for a single meson decaying into two nucleons). We can now divide out 
by T to get the transition function per unit time. But we still have to worry about 
summing over all final states. There are two steps: the first is to integrate over all 
possible momenta of the final particles: V J d^Pi/{2TT)^. The factors of spatial volume 
V in this measure cancel those in (3.87), while the factors of l/2Ep. in (3.87) conspire 
to produce the Lorentz invariant measure for 3-momentum integrals. The result is an 
expression for the density of final states given by the Lorentz invariant measure 

dn=(2.)V')(p.-p,) n ||4 

final states 

The second step is to sum over all final states with different numbers (and possibly 
types) of particles. This gives us our final expression for the decay probability per unit 
time, r = P. 

r = ^ E / l-^/^l'^n (3.89) 

final states 

r is called the width of the particle. It is equal to the reciprocal of the half-life r = l/F. 
3.6.3 Cross Sections 

Collide two beams of particles. Sometimes the particles will hit and bounce off each 
other; sometimes they will pass right through. The fraction of the time that they collide 
is called the cross section and is denoted by a. If the incoming fiux F is defined to 
be the number of incoming particles per area per unit time, then the total number of 
scattering events N per unit time is given by, 

N = Fa (3.90) 

We would like to calculate cr from quantum field theory. In fact, we can calculate a more 
sensitive quantity da known as the differential cross section which is the probability 
for a given scattering process to occur in the sohd angle {9, (f)). More precisely 

Differential Probability 1 1 , - ,9 

d(7 = L4fiPdn (3.91) 

Unit Time x Unit Flux AEiE^V F^^'^ ^ ' 

where we've used the expression for probability per unit time that we computed in the 
previous subsection. E\ and E2 are the energies of the incoming particles. We now 
need an expression for the unit flux. For simplicity, let's sit in the center of mass frame 
of the collision. We've been considering just a single particle per spatial volume V ^ 
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meaning that the flux is given in terms of the 3-velocities Vi as F = \vi — V2\/V. This 
then gives, 



da 



1 



1 



Afil'^dU 



(3.92) 



AE1E2 \v\ — V2 



If you want to write this in terms of momentum, then recall from your course on special 
relativity that the 3-velocities Vi are related to the momenta by -tT = p/my/l — — 



Equation (3.92) is our final expression relating the S-matrix to the differential cross 
section. You may now take your favorite scattering amplitude, and compute the proba- 
bility for particles to fly out at your favorite angles. This will involve doing the integral 
over the phase space of final states, with measure dll. Notice that different scattering 
amplitudes have different momentum dependence and will result in different angular 
dependence in scattering amplitudes. For example, in 0^ theory the amplitude for tree 
level scattering was simply A = —A. This results in isotropic scattering. In contrast, 
for nucleon-nucleon scattering we have schematically A ^ {t — m^)~^ + {u — m^)~^. 
This gives rise to angular dependence in the differential cross-section, which follows 
from the fact that, for example, t = — 2|p|^(l — cos^), where 6 is the angle between the 
incoming and outgoing particles. 

3.7 Green's Functions 

So far we've learnt to compute scattering amplitudes. These are nice and physical (well 
- they're directly related to cross-sections and decay rates which are physical) but there 
are many questions we want to ask in quantum field theory that aren't directly related 
to scattering experiments. For example, we might want to compute the viscosity of 
the quark gluon plasma, or the optical conductivity in a tentative model of strange 
metals, or figure out the non-Gaussianity of density perturbations arising in the CMB 
from novel models of inflation. All of these questions are answered in the framework of 
quantum fleld theory by computing elementary objects known as correlation functions. 
In this section we will briefly deflne correlation functions, explain how to compute them 
using Feynman diagrams, and then relate them back to scattering amplitudes. We'll 
leave the relationship to other physical phenomena to other courses. 

We'll denote the true vacuum of the interacting theory as We'll normalize H 
such that 



(3.93) 
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and (n| n) = 1. Note that this is different from the state we've called |0) which is the 
vacuum of the free theory and satisfies i^o |0) = 0. Define 

)^{n\T(PH{xi)...(f>HM\^) (3.94) 

where (f)H is (f) in the Heisenberg picture of the full theory, rather than the interaction 
picture that we've been dealing with so far. The G*^"^ are called correlation functions, 
or Green's functions. There are a number of different ways of looking at these objects 
which tie together nicely. Let's start by asking how to compute G^"^ using Feynman 
diagrams. We prove the following result 

Claim: We use the notation (pi = (/){xi), and write (pin to denote the field in the 
Heisenberg picture, and (pu to denote the field in the interaction picture. Then 

G^-\xi, . . . , x„) = (0| r0,H . . . \n) = (3.95) 

where the operators on the right-hand side are evaluated on |0), the vacuum of the free 
theory. 

Proof: Take ti > t2 > ■ ■ ■ > tn- Then we can drop the T and write the numera- 
tor of the right-hand side as 

(O|t/7( + OO,ti)0ut/(ti,t2)02/ ... 0„/t//(t„,-Oo)|O) 

= (0| Ui{+00, h)(j)iH ...(pnH Ui{tn, -Oo) |0) 

where we've used all the intermediate factors of Ui{tk,tk+i) — Texp{—i J^^^^ Hj) to 
convert (pj into (pH- Now let's deal with the two remaining C/(i,±oo) at either end of 
the string of operators. Consider an arbitrary state j^^) 

(^1 Ui{t, -oo) |0) = (^1 U{t, -oo) |0) (3.96) 

where U{t, — oo) is the Schrodinger evolution operator, and the equality above follows 
because i^o |0) = 0. Now insert a complete set of states, which we take to be energy 
eigenstates oi H — Hq + i^int, 



(*|C/(i,-oo)|0) = (*|C/(i,-oo) 



|0) + (n| 



|0) 



(^1 Jl) 0) + lim y e'^"(*'-*)(*|n) (n|0) (3.97) 



t'->-oo ' 
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But the last term vanishes. This follows from the Riemann-Lebesgue lemma which says 
that for any well-behaved function 




lim / dxf{x)e'^''' = 



(3.98) 



Why is this relevant? The point is that the in (3.97) is really an integral J dn, 
because all states are part of a continuum due to the momentum. (There is a caveat 
here: we want the vacuum |Q) to be special, so that it sits on its own, away from the 
continuum of the integral. This means that we must be working in a theory with a 
mass gap - i.e. with no massless particles). So the Riemann-Lebesgue lemma gives us 



(Notice that to derive this result, Peskin and Schroeder instead send t — > — oo in a 
slightly imaginary direction, which also does the job). We now apply the formula 
(3.99), to the top and bottom of the right-hand side of (3.95) to find 



which, using the normalization f2) = 1, gives us the left-hand side, completing the 



3.7.1 Connected Diagrams and Vacuum Bubbles 

We're getting closer to our goal of computing the Green's functions G^"^ since we 
can compute both (0| T0/(a;i) . . . ^/(a;^) 5 |0) and (Ol^lO) using the same methods 
we developed for S-matrix elements; namely Dyson's formula and Wick's theorem or, 
alternatively, Feynman diagrams. But what about dividing one by the other? What's 
that all about? In fact, it has a simple interpretation. For the following discussion, 
we will work in 0^ theory. Since there is no ambiguity in the different types of line 
in Feynman diagrams, we will represent the particles as solid lines, rather than the 
dashed lines that we used previously. Then we have the diagramatic expansion for 



lim {^\u{t,t') |o) = {n\o) 



(3.99) 



{0\Q) {Q\mH...CpnH\n) {n\0) 

(o|Q) {n\o) 



(3.100) 



proof. 



{0\S\0). 




(3.101) 



These diagrams are called vacuum bubbles. The combinatoric factors (as well as the 
symmetry factors) associated with each diagram are such that the whole series sums 



to an exponential, 



(0|5|0)=exp ^ g + O + + ...^ (3.102) 



So the amplitude for the vacuum of the free theory to evolve into itself is (0| 5 |0) = 
exp{all distinct vacuum bubbles). A similar combinatoric simplification occurs for generic 
correlation functions. Remarkably, the vacuum diagrams all add up to give the same 
exponential. With a httle thought one can show that 

(0| r0i . . . 0„ |0) = (^Y^ connected diagrams) (0| -5 |0) (3.103) 

where "connected" means that every part of the diagram is connected to at least one 
of the external legs. The upshot of all this is that dividing by (0| 5* |0) has a very nice 
interpretation in terms of Feynman diagrams: we need only consider the connected 
Feynman diagrams, and don't have to worry about the vacuum bubbles. Combining 
this with (3.95), we learn that the Green's functions G^"\xi . . . , x„) can be calculated 
by summing over all connected Feynman diagrams, 

{Q\T(I)h{xi) . . . (t)H{xn) \^) = ^ Connected Feynman Graphs (3.104) 

An Example: The Four-Point Correlator: (Vt \ T(f)H{xi) . . . (pnixi) 

As a simple example, let's look at the four-point correlation function in 0^ theory. The 
sum of connected Feynman diagrams is given by, 

Xi X2 

° ° + 2 Similar + + V ^ ^"^^^^ + ... 

X3° °X4 X3=^ ^X4 X3v/\x4 



All of these are connected diagrams, even though they don't look 
that connected! The point is that a connected diagram is defined 
by the requirement that every line is joined to an external leg. An 
example of a diagram that is not connected is shown in the figure. 
As we have seen, such diagrams are taken care of in shifting the 
vacuum from |0) to \fl). 

Feynman Rules 

The Feynman diagrams that we need to calculate for the Green's functions depend on 
Xi, . . . , Xn- This is rather different than the Feynman diagrams that we calculated for 




Figure 20: 
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the S-matrix elements, where we were working primarily with momentum eigenstates, 
and ended up integrating over all of space. However, it's rather simple to adapt the 
Feynman rules that we had earlier in momentum space to compute G^"'\xi . . . ,a;„). 
For (f)^ theory, we have 

• Draw n external points xi, . . . , Xn, connected by the usual propagators and ver- 
tices. Assign a spacetime position y to the end of each line. 

• For each line ^ y from x to y write down a factor of the Feynman propagator 

^f{x - y). 

• For each vertex position y, write down a factor of —iX J dt^y. 
3.7.2 Prom Green's Functions to S-Matrices 

Having described how to compute correlation functions using Feynman diagrams, let's 
now relate them back to the S-matrix elements that we already calculated. The first 
step is to perform the Fourier transform, 

G'(")(xi,...,x„) (3.105) 

These arc very closely related to the S-matrix elements that we've computed above. The 
difference is that the Feynman rules for G''-"-'(,Ti, . . . , effectively include propagators 
Ap for the external legs, as well as the internal legs. A related fact is that the 4- 
momcnta assigned to the external legs is arbitrary: they arc not on-shcll. Both of these 
problems are easily remedied to allow us to return to the S-matrix elements: we need 
to simply cancel off the propagators on the external legs, and place their momentum 
back on shell. We have 

n' n 

{p[, ...,p'JS-l\p,...,Pn)^ {-zr+^' Hip:' - m') n(p? - m') (3.106) 

i=i j=i 

X G'^-+-'\-p[,...,-p',„p,,...,Pn) 

Each of the factors {p^—irP') vanishes once the momenta are placed on-shcll. This means 
that we only get a non-zero answer for diagrams contributing to G^^\xi, . . . , Xn) which 
have propagators for each external leg. You might think they all do, but it's not true! 
Only diagrams that are fully connected, meaning each external point is connected to 
each other external point, have this property. For example, of the diagrams that we 
wrote down which contribute to the four-point function (Q| T(j)H{xi) . . . 4>h{xa) |^), only 
^^^^ will survive the multiplication by on-shell propagators in (3.106) to contribute to 
the S-matrix for meson scattering in 0^ theory. 



G'(")(pi,...,p,) = j 
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So what's the point of all of this? We've understood that ignoring the connected 
diagrams is related to shifting to the true vacuum But other than that, introducing 
the Green's functions seems like a lot of bother for little reward. The important point 
is that this provides a framework in which to deal with the true particle states in 
the interacting theory through renormalization. Indeed, the formula (3.106), suitably 
interpreted, remains true even in the interacting theory, taking into account the swarm 
of virtual particles surrounding asymptotic states. This is the correct way to consider 
scattering. In this context, (3.106) is known as the LSZ reduction formula. You will 
derive it properly next term. 
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